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ABSTRACT 

The unexpected dimness of Type la supernovae at redshifts z < 1 has over the past 

■ 7 years been seen as an indication that the expansion of the universe is accelerating. 
A new model cosmology, the "fractal bubble model", has been proposed by one of us, 
based on the idea that our observed universe resides in an underdense bubble remnant 
from a primordial epoch of cosmic inflation, together with a new solution for averaging 
in an inhomogeneous universe. Although there is no cosmic acceleration, it is claimed 
that the luminosity distance of type la supernovae data will nonetheless fit the new 
model, since it mimics a Milne universe at low redshifts. In this paper the hypothesis is 

j> ■ tested statistically against the available type la supernovae data by both chi-square and 

Bayesian methods. While the standard model with cosmological constant S7a = 1 — 
is favoured by a Bayesian analysis with wide priors, the comparison depends strongly 
on the priors chosen for the density parameter, £l m . The fractal bubble model gives 
| better agreement generally for fi m < 0.2. It also gives reasonably good fits for all the 

range, f2 m = 0.01-0.55, allowing the possibility of a viable cosmology with just baryonic 

■ matter, or alternatively with both baryonic matter and additional cold dark matter. 

Subject headings: Cosmology: theory - Cosmology: large-scale structure of universe 
'• Cosmological parameters — Cosmology: dark matter — Cosmology: early universe 

C3 



^ ' 1. Introduction 



For nearly a decade it has been assumed that the universe is presently undergoing a period 
of acceleration driven by an exotic form of dark energy. In a new model proposed by one of us 
(Wiltshire 2005a), it has been claimed that it is possible to fit type la supernovae (SNela) luminosity 
distances without exotic forms of dark energy, and without cosmic acceleration. The new model 
was developed on the basis of a suggestion of Kolb, Matarrese, Notari and Riotto (2005), that 
cosmological evolution must take into account super-horizon sized remnant perturbations from an 
epoch of primordial inflation. 

The new model cosmology of Wiltshire (2005a, 2005b) - henceforth the Fractal Bubble (FB) 
Model - shares the basic assumption of Kolb et al. (2005) that the observed universe resides within 
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an underdense bubble, S, in a bulk universe which has the geometry of a spatially flat Friedmann- 
Robertson-Walker (FRW) geometry on the largest of scales. In other respects, however, the FB 
model differs substantially from that of Kolb et al. (2005). In particular, the FB model is based on 
exact solutions of Einstein's equations, whereas Kolb et al. (2005) based their analysis on linearized 
perturbations, with a number of approximations. Furthermore, while Kolb et al. suggested that the 
deceleration parameter would behave as q(r) — > — 1 at late times, in the FB model q(r) — > 0. The 
FB model has the added advantage that it can make quantitative predictions for many cosmological 
quantities in the epoch of matter domination based on two parameters, the Hubble constant, H Q , 
and the density parameter, Q rn . 

The FB model makes two crucial physical assumptions. The first, as just mentioned, is that 
our observed universe resides in an underdense bubble, which is a natural outcome of primordial 
inflation. The second assumption is that on account of the bottom-up manner in which structure 
forms and spatial geometry evolves in an inhomogeneous universe with the particular density per- 
turbations that arise from primordial inflation, the local clocks of isotropic observers in average 
galaxies in bubble walls at the present epoch, are defined in a particular way. The small scale 
bound systems which correspond to typical stars in typical galaxies retain local geometries with 
an asymptotic time parameter frozen in to match the cosmic time of the average surfaces of ho- 
mogeneity in their past light cones at the epoch they broke away from the Hubble flow. That 
local average geometry was a spatially flat Einstein-de Sitter universe, even though the presently 
observable universe, S, is underdense when averaged on the larger spatial scales where structure 
forms later. These arguments are clarified at length by Wiltshire (2005b). 

The physical implication of this is that we must consider an alternative solution to the fitting 
problem in cosmology (Ellis and Stoeger 1987). The open FRW geometry 
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retains meaning as an average geometry for the whole universe at late epochs, but most closely 
approximates the local geometry only in voids, these being the average spatial positions on spatial 
hypersurfaces. However, in the matching of asymptotic geometries between those scales and the 
non-expanding scales of bounds stars and star clusters in average galaxies in bubble walls, it is 
assumed that a non-trivial lapse function j(t) = ^ enters, so that observers in such galaxies 
describe the same average geometry (1) in terms of a different conformal frame 

dS 2 = 7 2 ( T )ds 2 , (2) 

where 



ds 2 = -dr 2 + a 2 (r) 



+ r 2 (dfl 2 + sin 2 fld0 2 ) 



(3) 



l + r 2 

and a = 7 _1 a. The lapse function is uniquely determined, and with a careful re-calibration of 
cosmic clocks, a model cosmology can be constructed, which has the potential to fit luminosity 
distances without cosmic acceleration. 
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It is the aim of this Letter to test this claim against the available SNela data. We adopt the 
notation and conventions of Wiltshire (2005a,2005b). 



2. Observable quantities 



Starting from the standard definition of the redshift z, 

-i i -^observed 

1 + Z= , 

^emitted 

the luminosity distance for the cosmological model given in (Wiltshire 2005a) is 

c(l + z)(2 + fi* 



(4) 
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where H Q is the currently measured value of the Hubble constant, ^ is a conveniently defined 
parameter, related to the current matter density parameter, Q m , according to 
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and rj is given by 



cosh ?? = —- + 
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Fig. 1. — Difference of distance modulus from that of an empty (Milne) universe, Afi, versus 
redshift. The dashed line corresponds to the ACDM model with Vt m = 0.27, = 0.73; the 
solid lines correspond to the Fractal Bubble Model (FBM); and the dotted lines corresponds the 
open, matter only, universe (FRW). For the FBM and FRW models, in each case the upper line 
corresponds to Q. m = 0.05, the lower line to Q m = 0.27. (H Q = 62.7 km sec^ 1 Mpc^ 1 in all cases.) 
For Q m = 0.05 FB line is indistinguishable from the horizontal axis over the range shown. 
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Fig. 2. — Confidence limits for the estimation of parameters for the ACDM model, with £l\ = 

i — Q m . 

3. Data analysis 

To make contact with observation the luminosity distance (5) is related to the distance modulus 
\i by the standard formula 

\i = m - M = 51og 10 (d L ) + 25. (8) 

Let us define A/x = m mo( joi — ^cmpty, where m. ompty is the apparent magnitude for the O m = 
universe. A plot of A/j, versus z for the FB model as compared to the ACDM and open FRW 
models is given in Fig. 1. For Q m = 0.05 the distance modulus of the FB model is indistinguishable 
from that of a Milne universe over the range < z < 2; while even the Q m = 0.27 FB model is 
closer to a Milne universe than the Q. m = 0.05 open FRW model. 

We have compared the new FB model to the supernova data using the "Gold data set" recently 
published by Riess et al. (2004). The \ 2 an d Bayes factor method employed followed the one 
outlined in Ng and Wiltshire (2001). Use of the luminosity distance (5) requires a series expansion 
near f2 Q = to avoid numerical problems: for fi Q < 0.02 we used a Laurent series to order 9. The 
results are presented here in Figs. 2 and 3, where the la, 2a and 3a confidence contours are plotted 
in the (Q m ,H Q ) parameter space. 

A good fit for the ACDM model was at H Q = 65 km sec _1 Mpc _1 and Q m = 0.27 with 
(x 2 = 178) (Riess et al. 2004). For the Fractal Bubble Model the best fit values are at H Q = 
^7(2^ km sec' 1 Mpc" 1 and include the la range ft m < 0.44 { X 2 = 191.7). The value 
of x 2 is affected so little by the value of VL m in this range that we do not quote a single best-fit 
value. The 2a confidence limit is O m < 0.55. Q m can be made arbitrarily small within the la 
bounds. Since recalibrated primordial nucleosynthesis bounds for baryonic matter are typically in 
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the range ~ 0.02-0.13, we see that the supernova data fit well for the FB model regardless of 
whether or not there is cold dark matter in addition to baryonic matter. 

A Bayes factor for the two models can be calculated from the data by 

R = [ dH o I exp(-xLdoi i/ 2 ) (q) 
12 JdH Jdn m exp(- X l odcl2 /2)- U 

This result can be interpreted using Table 1. For very wide priors, namely 0.001 < Q m < 0.5 and 
55 < H Q < 75 km sec^ 1 Mpc" 1 , the Bayes factor is 396 in favour of ACDM. However, we note that 
the Bayes factor is very sensitive to the range chosen for the priors, since the best-fit values for the 
two models are include quite distinct values of Cl m . This is seen in Table 2, where narrower priors 
are used for H , and a variety of priors for Q, m . The examples with baryonic matter only use values 
of the baryon matter density from primordial nucleosynthesis bounds, as specifically recalibrated 
for the FB model by Wiltshire (2005b): it is found that the density fraction of ordinary baryonic 
matter is 2-3 times that which is conventionally obtained, and typically f2&~0.1 is expected. 



Table 1: Interpretation of Bayes Factors 



B\2 


Strength of evidence for H± over H2 


1 to 3 


Not worth more than a bare mention 


3 to 20 


Positive 


20 to 150 


Strong 


> 150 


Very Strong 



Table 2: Bayes Factor comparison of a Qa = 1 — ACDM model versus the FB model with a 
fixed prior for H = 100/ikm sec -1 Mpc -1 , 58 < h < 72, and varying priors for Q m . The baryonic 
density parameter bounds for the FB model are determined by Wiltshire (2005b). 



Q m prior 


Model type 


-B(ACDM):(FBM) 


Model favoured 


0.01 <n m < 0.5 


Wide priors 


396 


ACDM 


0.2 < n rn < 0.3 


favoured CDM range 


649 


ACDM 


0.3 < a rn < 0.5 


high density CDM 


1014 


ACDM 


0.01 < n m < 0.2 


low density CDM 


0.53 


FBM (slightly) 


0.08 < n m < 0.13 


baryonic matter, low D/H 


2.6 x 10" 5 


FBM 


0.02 < n m < 0.06 


baryonic matter, high D/H 


1.3 x 10~ 14 


FBM 



To characterise the parameter space better, we also give a plot of the x 2 value as a function of 
Q. m for the best-fit value of H Q in each model in Fig. 4, and similarly plot Xfbm/Xacdm f° r the two 
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models for the whole parameter space spanned by 0, m and H Q in Fig. 5. The Fractal Bubble Model 
is favoured for lower values of £l m and H , and ACDM for higher values of both parameters. Also 
observe from Fig. 4 how insensitive the FB model is to the value of Q. m in the range < fl m < 0.5. 



4. Data reduction 

While SNela serve as excellent standard candles there is some dispersion observed in their 
absolute magnitudes. However, a relation first proposed by Phillips (1993), relates the absolute 
magnitude, M, to the change in apparent magnitude from peak to 15 days after peak, Am 15 , as 
measured by clocks in the rest frame of the supernova. Various methods have been used to fit light 
curves to obtain absolute magnitudes (Drell, Loredo & Wasserman 2000, Tonry 2003). 

It is also well demonstrated that redshift causes a broadening of light curves (Leibundgut 1996) 
corresponding to the phenomenon of cosmic time dilation. In standard cosmology the redshift is 
entirely attributable to cosmic expansion. It is then a simple matter to transform the observed 
light curve to the rest frame. In the FB model 

1 + , = !%^, (io) 

a(T cm it)7 

so that the measured redshift differs from that inferred from (1). However, for null geodesies the 
transformation to the rest frame only involves only the overall (1 + z) factor of the measured 
redshift, so that the light curve broadening determination will be unaffected. Only measurements 
which involved some independent determination of 7(r cmit )/7 , would lead to measurable differences 
between this model and models based on the conventional identification of clock rates. 

It is impossible to comment on the Phillips relations, and other related aspects of data-fitting, 
without access to the raw data. However, since these are relations are empirical, it is not likely 
that there would be significant differences in the FB model. 

The only likely systematic errors that are as yet unaccounted for in reducing the data for 
comparison with the FB model are those effects directly related to the inhomogeneous geometry. 
The metric (1) represents an average geometry only, and for the redshifts which presently constitute 
the bulk of the sample, a better understanding of intermediate scales in the fitting problem may 
be required. Studies of local voids - at scales well within the bubble S~ do reveal that measurable 
effects on cosmological parameters are possible (Tomita 2001). In the present model, the average 
rate of expansion of typical 30/i _1 Mpc voids would differ from that of much larger voids. While the 
dominant effect on luminosity distances suggested here is a new effect due to the identification of an 
alternative homogeneous time parameter in solving the "fitting problem" , further scale dependent 
corrections in the fractal geometry are certainly to be expected. 
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5. Conclusion 

From the standard Bayesian analysis with wide priors we would conclude that on the basis of 
the "Gold data set" of (Riess et al. 2004) , the standard ACDM model is very strongly favoured over 
the FB model. Nonetheless, this conclusion depends strongly on the priors assumed, as illustrated 
by Table 2 and Figs. 4 and 5. We must also bear in mind that the average geometry (l)-(3) is just 
a first approximation to the fitting problem. A better understanding of the differential expansion 
rates of different characteristic sizes of voids in the inhomogeneous geometry may improve the fit 
of the FB model at those redshifts which are most significant for the claims of cosmic acceleration, 
rather than a nearly empty universe deceleration. The FB model has other significant advantages 
in terms of it greatly increased expansion age (Wiltshire 2005a) , which has the potential to explain 
the formation of structure at large redshifts. Thus it should be considered as a serious alternative 
to models with dark energy, quite apart from its intrinsic appeal of being a model based simply 
on general relativity and primordial inflation without needing exotic dark energy, the presence of 
which at the present epoch would be a mystery to physics. 

In the end, nature is the final arbiter, and future supernova observations at higher redshifts, 
such as those of the SNAP mission, will certainly be able to much more decisively distinguish 
between the FB model and ACDM, given the significantly more rapid deceleration of ACDM at 
redshifts z > 2. 

Acknowledgements: This work was supported by the Marsden Fund of the Royal Society of New 
Zealand. We thank Michael Albrow for helpful comments. 

Note added: After this work appeared an independent Snela data analysis of the Milne universe 
- considered as an approximation to a variety of alternative theories - has been performed (Sethi, 
Dev & Jain 2005). They obtain a comfortable fit just as we would expect from our results given 
that the FB model mimics a Milne universe for z < 2. 
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Fig. 3. — Confidence limits for the estimation of parameters for the Fractal Bubble Model . 
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Fig. 4. — A slice of the confidence plot for the Fractal Bubble Model at its best-ht value H Q = 62.683 
(solid line) for < Q m < 0.5, and for the ACDM model at its best-ht value H = 64.145 (dashed 
line). 
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Fig. 5. — A contour plot of Xfbm/Xacdm ^ or -Fractai Bubble Model vs the ACDM model. The 
Fractal Bubble Model fares better in the lower left region (white), while the ACDM model is 
favoured in the upper right region (lightly shaded). 



